2D stationary resistive MHD flows: borderline 
to magnetic reconnection solutions 



D.H. Nickeler a '*, H.-J. Fahr b 



a Astronomical Institute, Utrecht University, Princetonplein 5, 3584 CC Utrecht, 

the Netherlands 

h Institut fur Astrophysik und Extraterrestrische Forschung, Auf dem Hiigel 71, 

53121 Bonn, Germany 



Abstract 

We present the basic equations for stationary, incompressible resistive MHD flows 
in two dimensions. This leads to a system of differential equations for two flux 
functions, one elliptic partial differential equation (Grad-Shafranov-like) for the 
magnetic flux function and one for the stream function of the flow. In these equa- 
tions two potentials appear: one potential is a generalized pressure. The second 
potential couples the magnetic and the flow shear components of the system. With 
the restriction to flux or at least line conserving flows one has to solve a modified 
Ohm's law. For the two dimensional case these are two coupled differential equa- 
tions, which represent the borderline between the resistive but flux conserving (or 
line conserving) case, and that of reconnective solutions. We discuss some simplified 
solutions of these equations. 

Key words: magnetic reconnection, resistive MHD, heliopause 



1 Introduction 



A general problem of astrophysics, also reappearing in heliospheric physics, is 
the structure of the contact discontinuity which is located between the reverse 
shock 1 of a stellar wind and the counterstreaming interstellar medium. The 
contact discontinuity is called the heliopause or more generally the astropause. 
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Here not only two different flows are encountering, but also the electromagnetic 
fields transported by the flows. The aim of our work is to find magnetic field 
configurations in the vicinity of the stag nation point i n the framework of MHD. 
Similar work has already been done bv lPriest et a l. (1994). But their analysis 
used linearizations of the MHD-equations and their model was restricted to 
constant resistivities. 



2 Governing equations for stationary resistive MHD— flows in 2D 



Introducing the streaming vector w = ^fpv with V • v — and the electric 
potential <j) e (x,y) one can write the resistive MHD equations in the following 
form 
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R z are the components of the generalized resistivity with Ohm's law given 
by E + v x B = R, and p is the plama pressure. The potentials ( and a 
are functions of x and y only and are the streaming function respectively 
the magnetic flux function. The potential S ensures the vanishing of the z- 
component of the Lorentz force, which guarantees the vanishing of the z- 
component of the pressure gradient and therefore that the equilibrium is two- 
dimensional. Thus, the potential S induces the appropriate coupling between 
w z and B z . The equ ations for incompressible s tationary resistive MHD in 2D 
had been derived by Neukirch fc Priest (1996), but without the relations for 
the components of the magne tic field and streaming vector in the invariant 
directi on, which are derived in|Nickeler fc Our procedure follows 

mainly Neukirch fc Priest ( 19961 ) . however, we admit here of the ^-component 
of the magnetic field and streaming vector. 

To ensure the condition for magnetic field line freezing into the bulk plasma 
flow, one has to obey Ohm's law and the induction equation 
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Fig. 1. Streamlines of a standard stagnation point of a symmetric astropause. The 
astropause is the y-axis. 

Derivations of this relation can be found in Priest fc Forbes! <[200rt and refer- 
ences therein and in Hesse fc Schindlerl fll988h . For A = this reduces to the 
well-known frozen- in field condition (flux conserving flow). 

In the stationary two-dimensional case the above criterion must be written as 
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for the stationary two dimensional case. Equa tion (6) is equiva lent to the 
critereon for line-freezing, which is given, e.g. bv lVasvliuna^ (|l972h . The term 
A(a) is the difference between the resistive term and the constant electric field 
£0 7^ in ^-direction 2 . A simple case is A(a) = constant = A = — E Q which 
neglects resistive terms. Assuming a standard stagnation point given by the 
flow field C = axy, which is shown in Fig. 1, we are able to calculate a wide 
class of solutions of ideal Ohm's law for the magnetic field by using the theory 
of characteristics. These solutions of ideal Ohm's law are given by 



a = ~~~ \/~p (j^i \n(x/x f + (1 + h) \n(y/y ) 2 ) 
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Field lines of a special solution of the Ohm's law around a standard stagnation 
point can be seen in Fig. 2. where in this special case the constants a, xq and yo 



Eq = would imply that the flow is purely field aligned 



3 




Fig. 2. Magnetic field lines resulting from solving Ohm's law with respect to the 
velocity field of Fig. 1 for Eq ^ 0; for Eq = the magnetic field would be exactly 
flow-aligned, if the shear components are zero. 

are set to 1 and k\ is set to -1, which could represent the magnetic fields in the 
vicinity around a stagnation point of an astrosphere, where the flow vector w 
vanishes at the "nose" (i.e. stagnation point) of the astropause. However, in the 
case of ideality (A = const.) the velocity field diverges, although \w\ remains 
finite, so that the "stagnation point", is only a "topological" stagnation point. 
Only for the field w the origin is a stagnation point. Here we have to choose 
the density to be zero on the separatrices to avoid magnetic singularities. 
Therefore around the separatrices there must be a kind of a "vacuum gap" to 
ensure that the mass current density vector y/pw remains finite. This leads to a 
double separatrix structure of the magnetic field (see Fig. 2). So the constraint 
of "ideainess" leads to a magnetic structure which is more complex then the 
underlying structure of the flow field. 



3 Problems and Discussion 



The class of solutions we found does not ensure that these are also solutions 
of the Euler equation. It seems to be more probable, that stationary MHD 
equilibria, conserving magnetic flux, do not exist close to the standard stag- 
nation point, at least not without other forces, e.g. pressure anisotropies. It 
seems that to date only for a constant resistivity two dime nsional annihilation 



solut ions (without shear components) have been found ([Neukirch fc Priest 



1996). So our ansatz gives a possibility for further investigations, to find non 
flux conserving, but line conserving solutions A = A(a), so called annihilation 
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solutions, as magnetic flux is annihilated. 
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